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ABSTRACT 



Contagious parthenogenesis — a process involving rare functional males produced 
by a parthenogenetic lineage which mate with coexisting sexual females result- 
ing in fertile parthenogenetic offspring — is one of the most striking mechanisms 
responsible for the generation of new parthenogenetic lineages. Populations of the 
parthenogenetic diploid brine shrimp Artemia produce fully functional males in 
low proportions. The evolutionary role of these so-called Artemia rare males is, 
however, unknown. Here we investigate whether new parthenogenetic clones could 
be obtained in the laboratory by mating these rare males with sexual females. We 
assessed the survival and sex ratio of the hybrid ovoviviparous offspring from pre- 
vious crosses between rare males and females from all Asiatic sexual species, carried 
out cross-mating experiments between Fl hybrid individuals to assess their fertility, 
and estimated the viability and the reproductive mode of the resulting F2 offspring. 
Molecular analysis confirmed the parentage of hybrid parthenogenetic F2. Our study 
documents the first laboratory synthesis of new parthenogenetic lineages in Artemia 
and supports a model for the contagious spread of parthenogenesis. Our results 
suggest recessive inheritance but further experiments are required to confirm the 
likelihood of the contagious parthenogenesis model. 



Subjects Biodiversity, Evolutionary Studies, Genetics, Zoology 
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INTRODUCTION 

Parthenogenesis in animals has evolved through different molecular mechanisms that 
influence the initial genetic variability of parthenogenetic strains and therefore have 
important implications on their evolutionary success and persistence {Simon et at, 
2003). One of the most striking mechanisms responsible for the generation of new 
parthenogenetic lineages is contagious parthenogenesis (Simon et al, 2003; Schon, Martens 
& van Dijk, 2009). This involves a parthenogenetic lineage able to produce functional 
males, which mate with coexisting sexual females producing fertile parthenogenetic 
hybrid offspring. These new parthenogenetic lineages will combine genetic diversity of 
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the maternal sexual species and their paternal parthenogenetic ancestor, including the 
genetic fragments linked to the parthenogenesis (Simon etal, 2003; Tucker et al, 2013). 
This mechanism has been documented in aphids and parasitoid wasps (Schneider etal, 
2002; Sandrock & Vorburger, 2011; Delmotte et al, 2013), and most extensively in the 
Daphnia pulex species complex (Innes &Hebert, 1988; Poland, Colbourne & Lynch, 2005). 
In North American D. pulex parthenogenetic lineages, at least two distinct unrecombined 
haplotypes on chromosome VIII and IX are implied in the sex-limited meiosis suppression 
(Lynch et al, 2008; Lads et al, 2012; Tucker et al, 2013). These haplotypes leading to 
obligate parthenogenesis in D. pulex stem from a single recent event of hybridization 
with its sister taxon D. pulicaria (Xu et al, 2013; Tucker et al, 2013). Multiple new 
parthenogenetic lineages have arisen since this event as males produced by asexual lineages 
spread these parthenogenesis-inducing haplotypes by mating with sexual females. 

Artemia, an anostracan branchiopod commonly known as brine shrimp, is a typical 
inhabitant of hypersaline inland lakes and coastal lagoons and salterns. This genus includes 
sexual species and lineages of obligate parthenogenetic populations of diverse ploidy levels 
(Abatzopoulos, 2002), which makes it a good model system to investigate evolutionary 
transitions between reproductive systems. Parthenogenetic populations are restricted to 
the Old World where they co-occur with several sexual species in sympatry in various areas 
(Abatzopoulos, 2002; Agh et al, 2007; Abatzopoulos et al, 2009; Maccari et al, 2013). All 
strains of Artemia can reproduce either ovoviviparously, with the release of free-swimming 
nauplii broods when they complete their development in the ovisac (therefore, without 
a dormant phase), or oviparously with the production of broods of diapausing cysts 
(Browne, 1980; Abatzopoulos, 2002). 

In Artemia, both sexual and asexual females are heterogametic (ZW) (Stefani, 1963; 
Bowen, 1963; Bowen, 1965; De Vos et al, 2013). Diploid parthenogenetic lineages reproduce 
through automictic parthenogenesis, although the cytological details are controversial 
(Cuellar, 1987). It appears that diploidy restoration results in female offspring genetically 
identical to the mother barring mutation or recombination (Abreu-Grobois, 1987; Stefani, 
1960). Parthenogenetic diploid Artemia populations produce fully functional males in low 
proportions (Stefani, 1964; Bowen et al, 1978; MacDonald & Browne, 1987; Maccari et al, 
2013). Abreu-Grobois & Beardmore (2001 ) showed that rare males remain heterozygous 
at the same allozyme loci as their mothers, suggesting that rare males are produced as a 
result of rare ZW recombination events. These 'rare males' can generate viable offspring 
when crossed with females of sexual Asiatic species (Bowen et al, 1978; Cai, 1993; Maccari 
etal, 2013), to which they are closely related genetically (Munoz et al, 2010; Maniatsi et 
al, 2011; Maccari, Amat & Gomez, 2013), but they are reproductively isolated with other 
more distantly related species (MacDonald & Browne, 1987). However, the evolutionary 
role of rare males in the generation of Artemia parthenogenetic lineages is unknown 
(Maccari et al, 2013). The occurrence of contagious parthenogenesis has been suggested 
in light of the polyphyletic nature of maternal diploid parthenogenetic lineages (Maccari, 
Amat & Gomez, 2013), but we do not know if rare males are able to transmit partheno- 
genesis to their offspring, a requisite for contagious parthenogenesis. In an early study, 
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Bowen et al. (1978) crossed two parthenogenetic rare males, one from Yamaguchi (Japan) 
and the other one from Madras (India), with one sexual female of A. urmiana and one 
A. franciscana respectively, and concluded that parthenogenetic reproduction could not be 
transmitted through males because they failed to obtain parthenogenetic offspring either 
in hybrid Fl, F2 or F2 backcross. 

Laboratory generation and establishment of unisexual lineages can be a useful tool to 
complement phylogenetic approaches to identify the mechanism involved in the transition 
from sexual to parthenogenetic reproduction. However, most laboratory hybrids often 
exhibit low fertility and survival, or show deformation and abnormalities (Vrijenhoek, 
1989; Mantovani et al., 1996). In vertebrates, the first successful laboratory generation of a 
unisexual hybrid involved the origin of the hybridogenetic fish Poeciliopsis monacha-lucida 
through crosses of P. monacha females with P. lucida males (Schultz, 1973). Laboratory 
hybrids of hemiclonal European water frog R. esculenta (Rana ridibunda x Rana lessonae) 
show faster larval growth, earlier metamorphosis, and higher resistance to hypoxic 
conditions than their parental species and the equivalent hybrids in nature (Hotz et al, 
1999). More recently, Lutes et al. (201 1 ) generated self-sustaining tetraploid lineages of 
parthenogenetic lizards by pairing males of diploid sexual species Aspidoscelis inornata 
with females of the triploid parthenogenetic species Aspidocelis exsanguis. In invertebrates, 
the first laboratory generation of clonal hybrids in D. pulex was obtained by crossing males 
from obligately parthenogenetic clones with cyclically parthenogenetic females (Innes & 
Hebert, 1988). In addition, new lineages of thelytokous parthenogenetic lineages have been 
obtained in the wasp Lysiphlebus fabarum and in a South African honeybee, Apis mellifera 
capensis (Lattorff, Moritz & Fuchs, 2005; Sandrock & Vorburger, 201 1 ) . 

Here we assess the reproductive role of rare males and investigate whether new 
parthenogenetic clones could be produced in the laboratory as support for the contagious 
origin of parthenogenetic lineages in Artemia. For this purpose, ( 1 ) we assess the survival 
and sex ratio of the hybrid ovoviviparous offspring obtained from the previous crosses 
from Maccari et al. (2013 ) between rare males and four Asiatic sexual species, (2) we carry 
out cross-mating experiments between these Fl hybrid individuals to assess their fertility, 

(3) we estimate the viability and the reproductive mode of the resulting F2 offspring; 

(4) finally we demonstrate genetically that parthenogenetic F2 are indeed the descendants 
of the original crosses. This study shows that Artemia has the potential of generating 
parthenogenetic strains through contagious parthenogenesis. 

MATERIALS AND METHODS 
Populations and mating experiments 

In a previous study, we set up mating experiments between rare males from the 
diploid parthenogenetic Artemia population from Bagdad (Iraq, hereafter PD) and 
sexual females from Asiatic Artemia species to assess the fertility and the reproductive 
potential of rare males (Maccari et al., 2013). The females used were from the sexual 
Asiatic populations, A. urmiana from Koyashskoe Lake (Ukraine, URM), A. sinica from 
Yuncheng Lake (China, SIN), A. tibetiana from Lagkor Co Lake (Tibet, TIB) and Artemia 
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sp. from Kazakhstan (Artemia Reference Center code - ARC1039, unknown locality, 
KAZ). These interspecific crosses resulted in viable ovoviviparous and oviparous Fl 
offspring with similar or higher viability than controls (intraspecific sexual crosses) 
{Maccari et al., 2013). 

Survival rate, sex ratio and reproductive performance of hybrid 
generations 

For this study, live nauplii obtained from each ovoviviparous Fl hybrid brood were reared 
separately in jars containing brine at 80 gL~ 1 salinity, kept at 20-24 °C under mild aeration 
at a 12D:12L photoperiod and fed a mixture of Dunaliella sp and Tetraselmis sp. (1:1) 
microalgae every other day When animals showed signs of reproductive maturity they 
were counted and sexed to estimate survival rates (the proportion of F2 offpring per pair 
attaining adulthood) and sex ratio (the proportion of males in the F2 offspring per pair). 
For this procedure the animals were placed in Petri dishes with seawater and anaesthetized 
with a few drops of freshwater saturated with chloroform and examined carefully under 
a binocular microscope. We tested for deviations from a 50% sex ratio per cross and per 
pair using a Chi-square goodness of fit test (Pearson's statistic) (Wilson & Hardy, 2002). 
Statistical analyses were performed with SPSS v. 15.0 (SPSS Inc., Chicago, USA). 

Reproductive performance of the Fl hybrid individuals was evaluated in Fl x Fl cross 
fertility tests. For this purpose, 24 randomly size-matched hybrid Fl male-female pairs 
from each cross were transferred into separate small glass beakers (60 ml) under the culture 
conditions described above. Lifetime quantitative and qualitative reproductive outputs 
of each pair were monitored every other day during culture medium renewal events. 
For each paired Fl female we counted the number of unfertilized and fertilized broods, 
distinguishing the latter in oviparous and ovoviviparous broods. Eggs from unfertilised 
broods were identified as they are all smaller and white. In ovoviviparous offspring we also 
recorded the number of live and dead nauplii, and the number of abortive embryos (pale 
yellow-orange eggs). When oviparous offspring was produced, we counted the number 
of normally shelled diapausing cysts (pale grainy surface floating in 200 gl _1 brine), as 
opposed to abortive, abnormally shelled embryos (bright brown colour cysts sinking in 
200 gl" 1 brine) {Maccari et al, 2013). 

Emerged F2 hybrid nauplii were reared until maturity as described above. They were 
counted and sexed to estimate their survival rate and sex ratio in the F2 generation. Then, 
males and females were individually isolated in containers until their deaths to check if 
females could reproduce in isolation, as would be expected in parthenogenetic individuals. 
It is possible that some parthenogenetic females could be sterile; in this case, our procedure 
will underestimate the frequency of parthenogenesis. The proportion of parthenogenetic 
female offspring produced in each cross was tested against the expectations of 25% if 
governed by a recessive allele in a single gene using a Chi-square goodness of fit test. 
In addition, to test whether the different crosses produced the same percentage of 
parthenogenetic female offspring we used a Chi-square homogeneity test. 
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Paternity analysis of parthenogenetic F2 individuals 

(a) Microsatellite analysis 

The F2 hybrid generation resulting from crosses between rare males and sexual females 
from A. urmiana and Artemia sp. from Kazakhstan included parthenogenetic individuals. 
In order to rule out contamination and confirm that they were F2 individuals resulting 
from the original crosses, we screened three microsatellite loci, previously screened in the 
parental individuals in another study (Maccari etal., 2013), in the parthenogenetic F2 
animals obtained. Each microsatellite locus (Apdq02TAIL, Apdq03TAIL and Apd05TAIL) 
(Munoz et at, 2008) was amplified separately in PCRs performed as described in Maccari 
et al. (2013). Alleles were scored using the CEQ Fragment Analysis software (Beckman 
Coulter™) and checked manually. If F2 individuals had a paternal allele in any of the loci 
this would confirm that they were descendants of the diploid parthenogenetic rare males. 

(b) Maternal lineage 

The F2 resulting from the rare male x sexual female cross and Fl x Fl cross should 
carry the maternal DNA of the sexual strain. To establish the maternal lineage of the 
parthenogenetic F2 offspring, a 709-bp fragment of mitochondrial cytochrome c oxidase 
subunit I (COI) gene region was amplified in the parental (FO) individuals, in the Fl 
offspring and in the parthenogenetic F2 individuals. Total DNA was extracted and PCR was 
carried out as described previously (Maccari et al, 2013). PCR amplifications were sent to 
MACROGEN for sequencing, and the resulting electrophoregrams were checked by eye 
using CodonCode Aligner v. 3.5 (CodonCode Corporation, Dedham, MA). 

RESULTS 

Survival rate and sex ratio of F1 hybrid offspring 

A total of 102 ovoviviparous hybrid Fl broods produced by the crosses between each 
combination of sexual species with rare males (Maccari et al, 2013) were reared to 
maturity. The live nauplii obtained in each brood were morphologically normal. Survival 
rates to adulthood were over 50% in all Fl hybrid offspring (Fig. 1 ), and were highest in the 
Fl PD x SIN (80%), andlowestinFl PD x URMandFl PD x TIB (ca. 56%) (for the codes 
of the hybrid crosses see Fig. 1). The overall mean sex ratio of Fl offspring across pairs 
ranged from 49% males in Fl PD x KAZ cross to 53% males in Fl PD x TIB cross and did 
not significantly differ from 50% in any cross (Fig. 1). 

Reproductive performance of F1 hybrid offspring 

Prior to setting up the crosses, all females were isolated from males for two weeks to ensure 
that they could not reproduce in isolation (i.e., they were sexual females). No Fl females 
were able to reproduce when isolated from males. Then, a total of 24 mating pairs (Fl 
hybrid female x Fl hybrid male) were set up for each Fl produced in each combination of 
sexual species with rare males. As some individuals died before mating, the final number of 
experimental pairs ranged from 10 to 22 per cross, which produced a total of 173 fertile and 
92 infertile F2 hybrid broods (Table 1). Ovoviviparous and oviparous F2 offspring viability 
is shown in Fig. 2. The percentage of abortive embryos was high in all crosses (between 
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Figure 1 Survival rate and sex ratio (overall percentage of males) in the Fl hybrid offspring from 
Artemia rare males and Asiatic sexual females. Fl hybrids are from parental crosses between Anemia 
urmiana (URM), A. sinica (SIN), A. tibetiana (TIB), Artemia sp. from Kazakhstan (KAZ) and diploid 
parthenogenetic Artemia rare males (PD). Error bars are standard deviations. 



Table 1 Number of total, fertilized, ovoviviparous and oviparous broods in Fl Artemia hybrid off- 
spring. Fl hybrids are from parental crosses between Artemia urmiana (URM), Artemia sinica (SIN), 
Artemia tibetiana (TIB), Artemia sp. from Kazakhstan (KAZ) and diploid parthenogenetic Artemia rare 
males (PD). 



Cross 


Pairs 


Total 


Fertilized 


Ovoviviparous 


Oviparous 






broods 


broods 


broods 


broods 


Fl PD x KAZ 


18 


80 


42 


37 


5 


Fl PD x URM 


16 


48 


26 


22 


4 


Fl PD x TIB 


10 


33 


18 


4 


14 


Fl PD x SIN 


22 


104 


87 


40 


47 



70% and 90%), while the proportion of live nauplii in all hybrid ovoviviparous broods 
was low (from 5% to 25%). In oviparous broods, the proportion of properly shelled cysts 
ranged from 25% in F2 PD x TIB to 61% in F2 PD x URM. 

Survival rate and sex ratio of F2 hybrid offspring 

A total of 103 F2 ovoviviparous broods were recorded (Table 1), of which 35 broods 
from 27 pairs, characterized by the greatest number of nauplii, were followed to assess 
the survival rate and the sex ratio of the F2 offspring. F2 nauplii were morphologically 
normal but they had low survival rates when compared to Fl nauplii (Fig. 3). No F2 
hybrid offspring produced by the crosses between rare male and A. tibetiana survived to 
maturity. The F2 PD x KAZ had the highest survival rate, about 37%, followed by the F2 
PD x SIN (34%) and F2 PD x URM (24%). The overall mean sex ratio across pairs was 
significantly female-biased in F2 PD x KAZ and F2 PD x URM crosses (12% and 22% of 
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Figure 2 Reproductive traits (offspring quantity and quality) in F2 hybrids between Artemia rare 
males and Asiatic sexual females. The viability of ovoviviparous and oviparous broods is shown. Error 
bars are standard deviations. 




F2 PD X SIN 

(N=356) 



Figure 3 Survival rate and sex ratio (overall percentage of males) in the F2 hybrid offspring from 
Artemia rare males and Asiatic sexual females. F2 hybrids are from crosses between F 1 hybrid individu- 
als which are obtained in the crosses between Artemia urmiana (URM), A. sinica (SIN), A. tibetiana (TIB), 
Artemia sp. from Kazakhstan (KAZ) and diploid parthenogenetic Artemia rare males (PD). Error bars are 
standard deviations. Asterisks (P < 0.05) indicate significant differences from 50% sex ratio (Chi-square 
goodness of fit test was employed). 



males respectively; x 2 = 111.25 and/ 2 = 16.49, 1 df,p < 0.05), but was non-significantly 
different from 50% in the F2 PD x SIN (43% of males; / 2 = 0, 1 df, p < 0.05) (Fig. 3). 
Furthermore, we observed differences in the sex ratio of the F2 offspring among different 
pairs from the same cross, in particular for F2 PD x KAZ and F2 PD x URM crosses 
(see Table 2). In the cross F2 PD x KAZ, which higher sample sizes, one pair produced 
offspring with an even sex ratio (pair 3) while the remaining five pairs had were female 
biased offspring (see Table 2). 
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Table 2 Sex ratio and parthenogenetic females found in F2 PD x KAZ, F2 PD x URM and F2 PD 
x SIN Artemia offspring. Asterisks (P < 0.05) indicate significant differences from 50% sex ratio 
(number of males/total individuals) (Chi-square goodness of fit test was employed). All females obtained 
were isolated until their deaths to determine their mode of reproduction. 





Pair 


Females 


Males 


Total 


Sex ratio 

(%) 


Parthenogenetic 
females/analysed 
females 


Parthenogenetic 
females (%) 


F2 PD x KAZ 


1 


10 


0 


10 


0.00** 


3/10 


30 




2 


10 


2 


12 


16.67* 


1/10 


10 




3 


7 


8 


15 


53.33 


0/7 


0 




4 


20 


0 


20 


0.00** 


6/10 


60 




5 


68 


2 


70 


2.86** 


2/4 


50 




6 


31 


1 


32 


3.13** 


- 


- 


Total 




146 


13 


159 




12/41 


29.27 


F2 PD x URM 


1 


16 


3 


19 


15.79** 


0/16 


0 




2 


2 


4 


6 


66.67 


0/2 


0 




3 


2 


0 


2 


0.00 


0/2 


0 




4 


3 


1 


4 


25.00 


1/3 


33.33 




5 


2 


1 


3 


33.33 


- 


- 




6 


2 


0 


2 


0.00 


- 


- 




7 


13 


2 


15 


13.37** 


1/13 


7.69 


Total 




40 


11 


51 




2/36 


5.56 


F2 PD x SIN 


1 


15 


13 


28 


46.43 


0/15 


0 




2 


13 


24 


37 


64.86 


0/13 


0 




3 


6 


3 


9 


33.33 


0/6 


0 




4 


1 


3 


4 


75.00 


0/1 


0 




5 


14 


12 


26 


46.15 


0/14 


0 




6 


10 


10 


20 


50.00 


0/10 


0 




7 


20 


18 


38 


47.37 


0/20 


0 




8 


23 


24 


47 


51.06 


0/23 


0 




9 


30 


41 


71 


57.75 


0/30 


0 




10 


5 


8 


13 


61.54 


0/5 


0 




11 


16 


0 


16 


0.00** 


0/16 


0 




12 


7 


0 


7 


0.00** 


0/7 


0 




13 


4 


1 


5 


20.00 


0/4 


0 




14 


14 


21 


35 


60.00 


0/14 


0 


Total 




178 


178 


356 




0 


0 



Generation of hybrid parthenogenetic individuals 

Some females isolated from males of all F2 hybrid offspring analysed (when males were 
present) reproduced parthenogenetically in two of the three crosses. Specifically, 12 out 
of 41 isolated females (29.27%) were parthenogenetic in F2 PD x KAZ (four out of 
the five offspring analysed, Table 2), and two out of 36 (5.56%) isolated females in F2 
PD x URM (two of five offspring analysed, Table 2). The percentages of parthenogenetic 
female offspring in the F2 crosses were significantly different from each other (x 2 = 7.24, 
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1 df, p < 0.05), and only that one in F2 PD x KAZ did not differ significantly from the 
expectations of 25% (x 2 — 0-4, 1 df, p > 0.05) under expectations of a recessive allele in 
a single locus determining parthenogenesis. In all but one case, parthenogenetic females 
were produced in offspring with significantly female-biased sex ratios (Table 2). None of 
the 21 F2 PD x SIN offspring included females that could reproduce parthenogenetically. 

Paternity analysis 

In order to examine the parentage of newly generated hybrid parthenogenetic individuals 
we integrated the information from the mitochondrial COI and from microsatellites 
markers (Table 3 ) . Six of the 1 0 analysed females from pair 4 of the cross F2 PD x KAZ were 
parthenogenetic and produced F3 clones. As expected, all of them shared their mtDNA 
haplotype with their sexual grandmother, and amplified one paternal allele in the two 
informative microsatellite loci, confirming that they were the offspring of the rare male 
used in the crosses. The F3 generation was overall composed by females and by two rare 
males with the same genotype as their F2 mothers. 

The F2 offspring of two pairs from the crosses PD x URM (pairs 4 and 7), composed 
of three and 13 females respectively, included a parthenogenetic female that produced 
F3 parthenogenetic clones. In both cases, the F2 parthenogenetic female shared its COI 
haplotype with its sexual grandmother. In one cross, one paternal allele was detected 
in the F2 hybrid female at each of the three microsatellite loci; in the other cross, the 
parthenogenetic female inherited one paternal allele at the two informative loci. Most 
individuals of the F3 generation, composed of females and one rare male in both crosses, 
have the same genotype as their F2 mothers, with a few exceptions that lacked one of 
the maternal alleles, suggesting some level of recombination consistent with automixis 
parthenogenesis. 

DISCUSSION 

This study reports for the first time the laboratory generation of parthenogenetic Artemia 
lineages through hybridization via rare males, i.e., through contagious parthenogenesis 
(Simon et al, 2003), shedding light on the possible evolutionary role of parthenogenetically 
produced males and the genetic basis of parthenogenesis in this genus. 

Contagious parthenogenesis may have important evolutionary consequences as it 
results in the repeated generation of new asexual genotypes, increasing the genetic 
diversity in parthenogens. This may counteract the loss of asexual genotypes resulting 
from the accumulation of deleterious mutations (Muller's ratchet) or gene conversion 
(Tucker et al, 2013) and could contribute to the evolutionary success of parthenogenesis 
(Simon et al, 2003). 

The occurrence of contagious parthenogenesis relies on regular or occasional hy- 
bridization with absence of complete reproductive isolation between parthenogenetically 
produced males and closely related sexual females (Simon et al, 2003). In a previous 
study, we showed the absence of prezygotic isolation between rare males and Asiatic 
sexual Artemia species since these males often coexist in the same environment of a 
sexual species (Abatzopoulos et al, 2006; Agh et al, 2007; Agh et al, 2009; Shadrin, 
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Table 3 Mitochondrial cytochrome c oxidase subunit I (COI) and microsatellite loci analyses for 
parental individuals (FO) and for parthenogenetic F2 and F3 offspring obtained from the hybrid 
Artemia crosses. Genotypes for three microsatellite loci (allele sizes in base pairs) are shown. Diagnostic 
alleles, that is, alleles present in the rare male grandfather and not in the grandmother are shown in bold 
in the grandfather and in the F2 and F3 offspring. '0' indicates the presence of null alleles; 'm' indicates a 
rare male. COI haplotypes as named in GenBank are shown. KAZSEX03: GU591387; APD02: DQ426825; 
AUKOY02: KF707698; AUKOY01: KF707699. 





Sample code 


Apd02 


Apd03 


Apd05 


COI 


Rare male x 












Artemia sp. Kazakhstan 


FO (F-Kaz 8) 


233-233 


213-245 


0-0 


KAZSEX03 




FO (M-Iraq 8) 


233-242 


208-231 


115-0 


APD02 




F2-8-2-2 


233-233 


231-245 


0-0 


KAZSEX03 




F2-8-2-3 


233-242 


231-245 


0-0 


KAZSEX03 




F2-8-2-4 


233-242 


231-245 


0-0 


KAZSEX03 




F2-8-2-5 


233-242 


231-245 


0-0 


KAZSEX03 




F2-8-2-6 


242-242 


231-245 


0-0 


KAZSEX03 




F2-8-2-8 


233-242 


231-245 


0-0 


KAZSEX03 




F3-8-2-2-3 


233-233 


231-245 


0-0 


KAZSEX03 




F3-8-2-2-5 


233-233 


231-245 


0-0 


KAZSEX03 




F3-8-2-2-10 


233-233 


231-245 


0-0 


KAZSEX03 




F3-8-2-2-12m 


233-233 


231-245 


0-0 


KAZSEX03 




F3-8-2-6-3 


242-242 


231-245 


0-0 


KAZSEX03 




F3-8-2-6-4 


242-242 


231-245 


0-0 


KAZSEX03 




F3-8-2-6-5 


242-242 


231-245 


0-0 


KAZSEX03 




F3-8-2-6-7m 


242-242 


231-245 


0-0 


KAZSEX03 




F3-8-2-8-1 


233-242 


231-245 


0-0 


KAZSEX03 




F3-8-2-8-2 


233-242 


231-245 


0-0 


KAZSEX03 




F3-8-2-8-3 


233-242 


231-245 


0-0 


KAZSEX03 




F3-8-2-8-4 


233-242 


231-245 


0-0 


KAZSEX03 


Rare male x A. urmiana 


FO (F-Koy 15) 


233-281 


207-0 


170-0 


AUKOY02 




FO (M-Iraq 15) 


254-233 


216-231 


115-185 


APD02 




F2-15-8-A 


254-254 


207-216 


185 


AUKOY02 




F3-15-8-A-1 


254-254 


216 


185 


AUKOY02 




F3-15-8-A-4 


254-254 


207-216 


185 


AUKOY02 




F3-15-8-A-5 


254-254 


207-216 


185 


AUKOY02 




F3-15-8-A-6 


254-254 


207-216 


185 


AUKOY02 




F3-15-8-A-7m 


254-254 


207 


185 


AUKOY02 


Rare male x A. urmiana 


FO (F-Koy 16) 


248-0 


208-0 


90-90 


AUKOY01 




FO (M-Iraq 16) 


233-251 


216-230 


117-189 


APD02 




F2-16-7-4 


248-251 


0-0 


90-117 


AUKOY01 



{continued on next page) 
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Table 3 (continued) 



Sample code 


Apd02 


Apd03 


Apd05 


COI 


F3-16-7-4-1 


248-251 


0-0 


90-117 


AUKOY01 


F3-16-7-4-2 


248-251 


0-0 


90-90 


AUKOY01 


F3-16-7-4-3 


248-251 


0-0 


90-117 


AUKOY01 


F3-16-7-4-5 


248-251 


0-0 


90-117 


AUKOY01 


F3-16-7-4-7m 


248-251 


0-0 


90-117 


AUKOY01 



Anufriieva & Galagovets, 2012; Van Stappen et at, 2007; Van Stappen, 2008; Zheng & 
Sun, 2013), show normal pairing behaviour and are fully functional and capable of 
fertilizing eggs from females of sexual Asiatic Artemia species producing viable hybrid 
offspring (Maccari et at, 2013). Under laboratory conditions, each combination of 
sexual species with rare males produced morphologically normal, viable sexual hybrid 
Fl. Their survival rate to adulthood was over 50% for all the hybrid populations, a high 
value if compared to survival of Fl of intraspecific crosses of the different Artemia species 
(Browne & Wanigasekera, 2000). 

We found that females constitute approximately 50% of each Fl hybrid population, 
an even sex ratio that usually characterizes Artemia sexual populations, and this was 
confirmed by their inability to reproduce without males. These results ruled out a 
dominant gene as the genetic basis of parthenogenesis. Although all laboratory Fl lines 
were found to combine ovoviviparous and oviparous reproduction, we observed a strong 
reduction in the reproductive output in all crosses when compared with the reproductive 
performance of the parental crosses (Maccari et at, 2013). Ovoviviparous broods were 
mostly made up by abortive embryos (more than 80%) in all the crosses and live nauplii 
represented only 25% of the offspring in the cross F2 PD x SIN, and less than 10% in 
all the other crosses (F2 PD x KAZ, F2 PD x URM and F2 PD x TIB). Oviparity, the 
production of dormant encysted embryos that are resistant to extreme environmental 
conditions, was represented by a variable quantity of properly shelled embryos, only 
25% in the F2 PD x TIB increasing up to 61% in F2 PD x URM. Similarly, a decline in 
nauplii F2 production occurs in the interspecific crosses between A. tibetiana and A. sinica 
(Van Stappen etat, 2003). 

In contrast to the high survival rates of Fl hybrids, hybrid breakdown was evident 
in the F2 generation. Nauplii from the F2 generations had low survival rates and were 
completely inviable in the F2 PD x TIB generation. The lower fertility level of Fl 
laboratory populations and the reduced viability of F2 hybrid individuals suggest partial 
genetic incompatibility between parthenogenetic males and sexual females. However, 
the production of some viable offspring both in Fl and F2 in all hybrid crosses is not 
so surprising given the recent evolutionary origin of diploid parthenogenetic lineages 
(Holocene) (Muhoz et at, 2010; Maccari, Amat & Gomez, 2013). 

In two of the three F2 generations (F2 PD x KAZ and F2 PD x URM) we identified 
14 hybrid females that upon reaching maturity were capable of parthenogenetic 
reproduction. Surprisingly, these parthenogenetic females were produced by pairs yielding 
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strongly female biased F2 offspring. Genetic analysis confirmed the parentage of the 
parthenogenetic lineages found as the F2 individuals inherited the COI haplotype from 
the sexual grandmother but included some paternal alleles at nuclear markers, showing 
that they were the offspring of the rare male used in the crosses. Our results contrast with 
previous observations suggesting that rare males in the genus Artemia are not capable to 
transmit parthenogenesis-inducing alleles {Bowen et al., 1978). 

The production of parthenogenetic individuals only in the second generation, suggests 
that the parthenogenesis-inducing alleles are recessive in Artemia. A single-locus recessive 
inheritance of obligate parthenogenesis also occurs in Apis mellifem capensis and in 
Lysiphlebus fabarum (Sandrock & Vorburger, 2011; Lattorff, Moritz & Fuchs, 2005; Lattorff 
et al, 2007). This is in contrast with D. pulex, where the sex- limited meiosis suppression 
genes are dominant and the asexual clones arise in the first generation (Innes & Hebert, 
1988). If a single recessive locus was responsible for parthenogenesis and there was no 
differential viability in Artemia, a 25% of parthenogenetic females would be expected in the 
F2 generation. The proportion of isolated females that reproduced parthenogenetically 
differed between the crosses. In the cross F2 PD x KAZ, the overall proportion of 
parthenogenetic F2 females was 29.27%, not significantly different from 25%, whereas 
in the cross F2 PD x URM this was much lower (5.56%) and significantly different 
from the expectations for a single recessive locus. These results suggest either differences 
in the mechanism underlying parthenogenesis between populations, or increased 
incompatibilities between PD and URM resulting in viability differences linked to the 
putative locus associated to parthenogenesis. The latter is supported by the lower viability 
of F2 PD x URM nauplii. The finding of parthenogenetic females only in sex-biased 
broods suggests that the inheritance of parthenogenesis has a more complex genetic basis, 
however. Given that females are heterogametic (WZ) {Bowen, 1963; Bowen, 1965; Stefani, 
1963) and that Fl females are sexual, we can rule out complete sex-linkage (Z-linkage) of 
the parthenogenesis determining gene, otherwise parthenogenesis should be apparent in 
the Fl, given that all Fl females are WZ with their Z chromosome presumably inherited 
from their asexual father. Sex-biased sex ratios are not uncommon in hybrid offspring and 
can be due to the evolution of sex-ratio distorters and counter evolution of suppressor 
genes in different lineages (Hurst & Pomiankowski, 1991). Our data suggests an interaction 
between a sex ratio distorter (possibly sex-linked) and a parthenogenetically determining 
factor. Alternatively, the same gene determining parthenogenesis could act as a sex ratio 
distorter in heterozygous Fl females, increasing the likelihood of transmission of the W 
chromosome. Our results do not support differential male mortality, as there was no 
correlation between brood survival and sex ratio (data not shown). These interpretations 
must be taken with caution given the limitations of our experimental design and data, as 
we analysed F2 broods where there was a larger number of nauplii, the survival of the F2 
was low, and we cannot rule out some effect of differential sterility. These factors might 
have biased our conclusions regarding the genetic basis of parthenogenesis. Therefore, to 
fully understand the genetic basis of parthenogenesis in Artemia additional crosses and a 
large set of marker loci will be necessary. 
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The ability of sexual females of A. urmiana and Artemia sp. from Kazakhstan to generate 
parthenogenetic clones when crossed with rare males is not surprising, as the two main 
mitochondrial haplogroups of diploid parthenogenetic Artemia lineages are related to 
these species (Munoz et ah, 2010; Maniatsi et ah, 2011; Maccari, Amat & Gomez, 2013). 
However, the more distantly related A. sinica (Baxevanis, Kappas & Abatzopoulos, 2006; 
Hon et ah, 2006) did not produce any parthenogenetic offspring, despite high survival 
rate in the F2, suggesting that the specific genomic background affect the expression of 
the gene inducing parthenogenesis. Although repeated gene flow between sexual females 
and asexual males through contagious parthenogenesis would be expected to result in a 
regular emergence of asexual strains with diverse maternal origins, the fact that just two, 
possibly three, maternal origins of parthenogenetic lineages have been identified {Munoz et 
ah, 2010; Maniatsi et ah, 2011; Maccari, Amat & Gomez, 2013) indicate that the incidence 
of contagious parthenogenesis, if this is the mechanism of origin, must be extremely low in 
natural environments. Indeed, the rare males must be present in the population at the same 
time as the sexual females of the related species, and given that both parthenogenetic and 
sexual species often have different ecological requirements, they may overlap just during 
part of each season (Amat et ah, 1991; Ghomari et ah, 2011). In addition, the percentage 
of rare male production by diploid parthenogenetic females is very low, about 1-16 in 
1000 (Maccari et ah, 2013). Then, as the parthenogenesis occurs in the second generation 
(i.e., is based on a recessive trait), a Fl x Fl mating must occur for parthenogenesis to 
appear in the offspring. Finally, F2 survival is very reduced, overall making the origin of a 
parthenogenetic lineage an unlikely event in the wild. 

Our study is the first to generate new parthenogenetic lineages in Artemia by mating 
rare males from parthenogenetic genotypes with sexual females, providing evidence that 
contagious parthenogenesis can potentially occur in the genus Artemia. This conclusion 
does not rule out that other mechanisms (spontaneous origin or hybridisation) might 
have been also responsible for the origin of parthenogenetic lineages. Demonstration of 
contagious parthenogenesis as the mechanism underlying parthenogenesis in Artemia in 
the wild will necessitate the use of genomic tools. Further studies on hybrid fitness would 
be necessary to estimate the strength of reproductive isolation and to compare the repro- 
ductive performance of laboratory-produced parthenogenetic clones with the parental 
parthenogenetic strains. The origin of independently reproducing parthenogenetic clones 
in the laboratory raises the question of the survival of these clones when competing with 
sympatric sexual species. 

Given that many parthenogenetic organisms produce males occasionally (van der 
Kooi & Schwander, 2014) and such males are still able to maintain their functionality, 
the occurrence of contagious parthenogenesis could be more widespread than currently 
acknowledged. 

ACKNOWLEDGEMENTS 

We wish to thank Paul Nichols for his help with microsatellite screening and Monica 
Barbosa, Eva Becerro and Diana Guinot for their help with the laboratory experiments. 



Maccari et al. (2014), PeerJ, DO1 1 0.771 7/peerj.439 



■ 13/18 



PeerJ 



We thank Maria Jose Carmona for her constructive suggestions on a previous version 
of this manuscript. We thank the editor Tanja Schwander, and David Innes and two 
anonymous reviewers for their constructive comments that substantially improved the 
manuscript. 



ADDITIONAL INFORMATION AND DECLARATIONS 



Funding 

This study was funded by the Plan Nacional CGL2008-03277 project to FA, sponsored 
by the Spanish Government MICIN. AG was supported by a National Environment 
Research Council (NERC) Advanced Fellowship (NE/B50 1298/1). MM was supported 
by a fellowship of the JAE Program from CSIC and European Social Fund. The funders had 
no role in study design, data collection and analysis, decision to publish, or preparation of 
the manuscript. 

Grant Disclosures 

The following grant information was disclosed by the authors: 
MICIN Plan Nacional: CGL2008-03277. 

National Environment Research Council (NERC) Advanced Fellowship: NE/B50 1298/1. 
JAE-CSIC. 

European Social Fund. 

Competing Interests 

The authors declare there are no competing interests. 

Author Contributions 

• Marta Maccari conceived and designed the experiments, performed the experiments, 
analyzed the data, wrote the paper, prepared figures and/or tables, reviewed drafts of the 
paper. 

• Francisco Amat conceived and designed the experiments, contributed 
reagents/materials/ analysis tools, reviewed drafts of the paper. 

• Francisco Hontoria contributed reagents/materials/analysis tools, reviewed drafts of the 
paper. 

• Africa Gomez conceived and designed the experiments, contributed 
reagents/materials/ analysis tools, wrote the paper, reviewed drafts of the paper. 

Supplemental Information 

Supplemental information for this article can be found online at http://dx.doi.org/ 
10.7717/peerj.439. 

REFERENCES 

Abatzopoulos TJ. 2002. Artemia: basic and applied biology. Dordretch: Kluwer Academic 
Publishers. 



Maccari et al. (2014), PeerJ, DO1 1 0.771 7/peerj.439 



14/18 



PeerJ 



Abatzopoulos TJ, Agh N, Vanstappen G, Razavi-Rouhani SM, Sorgeloos P. 2006. Artemia sites in 

Iran. Journal of Marine Biology 86:299-307 DOI 10.1017/S0025315406013154. 
Abatzopoulos TJ, Amat F, Baxevanis AD, Belmonte G, Hontoria F, Maniatsi S, Moscatello S, 

Mura G, Shadrin NV. 2009. Updating geographic distribution of artemia urmiana Gilnther, 

1890 (branchiopoda: anostraca) in Europe: an integrated and interdisciplinary approach. 

International Review of Hydrobiology 94:560-579 DOI 10.1002/iroh.200911147. 
Abreu-Grobois FA, Beardmore JA. 2001. The generation of males by diploid parthenogenetic 

Artemia cannot occur in the way Stefani suggested. In: Proceedings of the 4th international large 

branchiopod symposium. La Paz, Baja California, Mexico, 1. 
Abreu-Grobois FA. 1987. A review of the genetics of Artemia. In: Sorgeloos P, Bengtson DA, 

Decleir W, Jaspers E, eds. Artemia research and its applications, vol. 1: morphology, genetics, 

strain characterization, toxicology. Wetteren: Universa Press. 
Agh N, Bossier P, Abatzopoulos TJ, Beardmore JA, Vanstappen G, Mohammadyari A, 

Ramihian H, Sorgeloos P. 2009. Morphometric and preliminary genetic characteristics 

of Artemia populations from Iran. International Review of Hydrobiology 94(2):194-207 

DOI 10.1002/iroh.200811077. 
Agh N, Abatzopoulos TJ, Kappas I, van Stappen G, Razavi Rouhani SM, Sorgeloos P. 2007. 

Coexistence of sexual and parthenogenetic Artemia Populations in Lake Urmia and 

neighbouring lagoons. International Review of Hydrobiology 92:48-60 

DOI 10.1002/iroh.200610909. 
Amat F, Hontoria F, Navarro JC, Gozalbo A, Varo 1. 1991. Bioecologia de Artemia (Crustacea, 

Branchiopoda) en la laguna de La Mata. Coleccion Tecnica Publ. Inst. Cultura "Juan Gil Albert". 

Diputacion de Alicante, p. 176. 
Baxevanis AD, Kappas I, Abatzopoulos TJ. 2006. Molecular phylogenetics and asexuality in the 

brine shrimp Artemia. Molecular Phylogenetics and Evolution 40:724-738 

DOI 10.1016/j.ympev.2006.04.010. 
Bowen ST, Durkin JP, Sterling G, Clark LS. 1978. Artemia hemoglobins: genetic variation in 

parthenogenetic and zygogenetic populations. The Biological Bulletin 155:273-287 

DOI 10.2307/1540952. 
Bowen ST. 1965. The genetics of Artemia salina V. Crossing over between the x and y 

chromosomes. Genetics 52:695-710. 
Bowen ST. 1963. The genetics of Artemia salina. II. White, a sex-linked mutation. Biological 

Bulletin \2A.:\1-Ti DOI 10.2307/1539564. 
Browne RA. 1980. Competition experiments between parthenogenetic and sexual strains of the 

brine shrimp, Artemia salina. Ecology 61:471-474 DOI 10.2307/1937409. 
Browne RA, Wanigasekera G. 2000. Combined effects of salinity and temperature on survival 

and reproduction of five species of Artemia. Journal of Experimental Marine Biology and Ecology 

244:29-44 DOI 10.1016/S0022-0981(99)00125-2. 
Cai Y. 1993. "OPI" Artemia and its crossing experiments. Journal of Ocean University of Qingdao 

23:52-58. 

Cuellar 0. 1987. The evolution of parthenogenesis: a historical perspective. In: Moens P, ed. 

Meiosis. Orlando, FL: Elsevier, 43-104. 
Delmotte F, Sabater-munoz B, Prunier-leterme N, Sunnucks P, Rispe C, Simon J, Url S. 2013. 

Phylogenetic evidence for hybrid origins of asexual lineages in an aphid species. Evolution 
57:1291-1303 DOI 10.1111/j.0014-3820.2003.tb00337.x. 



Maccari et al. (2014), PeerJ, DO1 1 0.771 7/peerj.439 



15/18 



PeerJ 



De Vos S, Bossier P, Van Stappen G, Vercauteren I, Sorgeloos P, Vuylsteke M. 2013. A first 

AFLP-based genetic linkage map for brine shrimp Artemia franciscana and its application in 

mapping the sex locus. PLoS ONE 8(3):e57585 DOI 10.1371/journal.pone.0057585. 
Eads BD, Tsuchiya D, Andrews J, Lynch M, Zolan ME. 2012. The spread of a transposon insertion 

in Rec8 is associated with obligate asexuality in Daphnia. Proceedings of the National Academy 

of Sciences of the United States of America 109(3):l-6 DOI 10.1073/pnas.lll9667109. 
Ghomari M, Selselet G, Hontoria F, Amat F. 2011. Artemia biodiversity in Algerian Sebkhas. 

Crustaceana 84(9): 1025-1039 DOI 10.1163/001121611X586729. 
Hotz H, Semlitsch RD, Gutmann E, Guex GD, Beerli P. 1999. Spontaneous heterosis in larval 

life-history traits of hemiclonal frog hybrids. Proceedings of the National Academy of Sciences of 

the United States of America 96:2171-2176 DOI 10.1073/pnas.96.5.2171. 
Hou L, Bi X, Zou X, He C, Yang L, Qu R, Liu Z. 2006. Molecular systematics of bisexual Artemia 

populations. Aquaculture Research 37:671-680 DOI 10.1111/j.l365-2109.2006.01480.x.. 
Hurst LD, Pomiankowski A. 1991. Causes of sex ratio bias may account for unisexual sterility in 

hybrids: a new explanation of Haldane's rule and related phenomena. Genetics 128:841-858. 
Innes DJ, Hebert PDN. 1988. The origin and genetic basis of obligate parthenogenesis in Daphnia 

pulex. Evolution 42:1024-1035 DOI 10.2307/2408918. 
Lattorff HMG, Moritz RFA, Crewe RM, Solignac M. 2007. Control of reproductive dominance by 

the thelytoky gene in honeybees. Biology Letters 3:292-295 DOI 10.1098/rsbl.2007.0083. 
Lattorff HMG, Moritz RFa, Fuchs S. 2005. A single locus determines thelytokous parthenogenesis 

of laying honeybee workers (Apis mellifera capensis). Heredity 94:533-537 

DOI 10.1038/sj.hdy.6800654. 
Lutes AA, Baumann DP, Neaves WB, Baumann P. 201 1. Laboratory synthesis of an 

independently reproducing vertebrate species. Proceedings of the National Academy of Sciences of 

the United States of America 108:9910-9915 DOI 10.1073/pnas.ll02811108. 
Lynch M, Seyfert A, Eads BD, Williams E. 2008. Localization of the genetic determinants of 

meiosis suppression in Daphnia pulex. Genetics 180:317-327 DOI 10. 1534/genetics. 107.084657. 
Maccari M, Amat F, Gomez A. 2013. Origin and genetic diversity of diploid parthenogenetic 

Artemia in Eurasia. PLoS ONE 8:e83348 DOI 10.1371/journal.pone.0083348. 
Maccari M, Gomez A, Hontoria F, Amat F. 2013. Functional rare males in diploid 

parthenogenetic Artemia. Journal of Evolutionary Biology 26:1934-1948 DOI 10.1 1 1 l/jeb.12191. 
MacDonald GH, Browne RA. 1987. Inheritance and reproductive role of rare males in a 

parthenogenetic population of the brine shrimp, Artemia parthenogenetica. Genetica 75:47-53 

DOI 10.1007/BF00056032. 
Maniatsi S, Baxevanis AD, Kappas I, Deligiannidis P, Triantafyllidis A, Papakostas S, 

Bougiouklis D, Abatzopoulos TJ. 2011. Is polyploidy a persevering accident or an adaptive 

evolutionary pattern? The case of the brine shrimp Artemia. Molecular Phylogenetics and 

Evolution 58:353-364 DOI 10.1016/j.ympev.2010.11.029. 
Mantovani B, Tinti F, Barilani M, Scali V. 1996. Current reproductive isolation between ancestors 

of natural hybrids in Bacillus stick insects (Insecta: Phasmatodea). Heredity 77:261-268 

DOI 10.1038/hdyl996.140. 
Murioz J, Gomez A, Green AJ, Figuerola J, Amat F, Rico C. 2010. Evolutionary origin and 

phylogeography of the diploid obligate parthenogen Artemia parthenogenetica (Branchiopoda: 

Anostraca). PLoS ONE 5:ell932 DOI 10.1371/journal.pone.0011932. 



Maccari et al. (2014), PeerJ, DO1 1 0.771 7/peerj.439 



■^16/18 



PeerJ 



Munoz J, Green AJ, Figuerola J, Amat F, Rico C. 2008. Characterization of polymorphic 

micro satellite markers in the brine shrimp Artemia (Branchiopoda, Anostraca). Molecular 

Ecology Resources 9:547-550 DOI 10.1111/j.l755-0998.2008.02360.x. 
Paland S, Colbourne JK, Lynch M. 2005. Evolutionary history of contagious asexuality in Daphnia 

pulex. Evolution 59:800-813 DOI 10.1111/j.0014-3820.2005.tb01754.x. 
Sandrock C, Vorburger C. 2011. Single-locus recessive inheritance of asexual reproduction in a 

parasitoid wasp. Current Biology 21:433-437 DOI 10.1016/j.cub.2011.01.070. 
Schneider MV, Beukeboom LW, Driessen G, Lapchin L, Bernstein C, Van Alphen JJM. 

2002. Geographical distribution and genetic relatedness of sympatrical thelytokous and 

arrhenotokous populations of the parasitoid Venturia canescens (Hymenoptera). Journal of 

Evolutionary Biology 15:191-200 DOI 10.1046/j.l420-9101.2002.00394.x. 
SchultzR. 1973. Unisexual fish: laboratory synthesis of a "species". Science 179:180-181 

DOI 1 0. 1 126/science. 1 79.4069. 1 80. 
Schon I, Martens K, van Dijk P. 2009. Schon I, Martens K, van Dijk P, eds. Lost sex: the 

evolutionary biology of parthenogenesis. Dordrecht: Springer. 
Shadrin NV, Anufriieva E, Galagovets E. 2012. Distribution and historical biogeography of 

Artemia Leach, 1819 (Crustacea, Anostraca) in Ukraine. International Journal of Artemia Biology 

2(2):30-42. 

Simon J-C, Delmotte F, Rispe C, Crease T. 2003. Phylogenetic relationships between 

parthenogens and their sexual relatives: the possible routes to parthenogenesis in animals. 
Biological Journal of the Linnean Society 79:151-163 DOI 10.1046/j.l095-8312.2003.00175.x. 

Stefani R. 1964. The origin of males in parthenogenetic populations of Artemia salina. Rivista di 
Biologia 57:147-156. 

Stefani R. 1963. La digametia femminile in Artemia salina Leach e la constituzione del corredo 

cromosomico nei biotipi diploide anfigonico e diploide partenogenetico. Caryologia 16:629-636 

DOI 10.1080/00087114.1963.10796102. 
Stefani R. 1960. LArtemia salina pathenogenetica a Cagliari. Rivista di Biologia 53:463-490. 
Tucker AE, Ackerman MS, Eads BD, Xu S, Lynch M. 2013. Population-genomic insights 

into the evolutionary origin and fate of obligately asexual Daphnia pulex. Proceedings 

of the National Academy of Sciences of the United States of America 110:15740-15745 

DOI 10.1073/pnas.l313388110. 
van der Kooi CJ, Schwander T. 2014. On the fate of sexual traits under asexuality. Biological 

Reviews 1 DOI 10.1111/brv.l2078. 
Van Stappen G, Sui L, Xin N, Sorgeloos P. 2003. Characterisation of high-altitude Artemia 

populations from the Qinghai-Tibet Plateau, PR China. Hydrobiologia 500:179-192 

DOI 10.1023/A:1024658604530. 
Van Stappen G, Yu H, Wang X, Hoffman S, Cooreman K, Bossier P, Sorgeloos P. 2007. 

Occurrence of allochthonous Artemia species in the Bohai Bay area, PR China, as confirmed 

by RFLP analysis and laboratory culture tests. Archiv fur Hydrobiologie 170(l):21-28 

DOI 10.1127/1863-9135/2007/0170-0021. 
Van Stappen G. 2008. Artemia biodiversity in Central and Eastern Asia. Ph.D. Thesis, Ghent 

University, Belgium, 224. 
VrijenhoekR. 1989. Genetic and ecological constraints on the origins and establishment of 

unisexual vertebrates. In: Dawley RM, Bogart JP, eds. Evolution and ecology of unisexual 

vertebrates. Albany, New York: New York State Museum, 24-31. 



Maccari et al. (2014), PeerJ, DO1 1 0.771 7/peerj.439 



H 17/18 



PeerJ 



Wilson K, Hardy ICW. 2002. Statistical analysis of sex ratios. In: Hardy I, ed. Sex ratios: concepts 

and research methods. Cambridge: University Press. 
Xu S, Innes DJ, Lynch M, Cristescu ME. 2013. The role of hybridization in the origin and spread 

of asexuality in Daphnia. Molecular Ecology 22:4549-4561 DOI 10.1111/mec.l2407. 
Zheng B, Sun SC. 2013. Review of the biogeography oiArtemia Leach, 1819 (Crustacea, Anostraca) 

in China. International Journal of Artemia Biology 3(l):20-50. 



Maccari et al. (2014), PeerJ, DO1 1 0.771 7/peerj.439 



18/18 



